Evolution of the crystal-field splittings in the compounds CeX (X=P, As, Sb, Bi), 
CeY (Y=S, Se, Te) and their alloys CeXi_^Y^ 
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The crystal-field splittings of the monopnictides and monochalcogenides of Cerium (CeX and 
CeY) and their alloys (CeXi-ajY^) are calculated by means of an ab initio many-body combined 
technique. The hybridization functions of the 4f states of Cerium with the conduction band for each 
material are obtained from first principles within the local density approximation (LDA) and are 
used as input for the Anderson impurity model, which is solved within a multi-orbital Non-Crossing 
Approximation (NCA). This realistic theoretical approach (LDA-NCA) is able to reproduce the 
experimental results for the crystal-field splittings of the CeX and CeY series in agreement with 
previous theoretical calculations. It is also able to describe the non-linear evolution of the splittings 
in the CeXi-^^Ya; alloys as a function of x. An analysis of the values of the crystal-field splittings 
in all the compounds can be done in depth in this contribution, due to a detailed knowledge of the 
band structure and crystal environment in combination with many-body physics. 



I. INTRODUCTION 

One of the important issues of the last years has been 
to produce realistic theoretical descriptions of strongly 
correlated electronic systems. These descriptions make 
a point out of including detailed electronic structure in- 
formation as well as many body dynamical effects on an 
equal footing. Examples of materials for which this kind 
of treatments becomes essential are Cerium compounds. 
In these systems, the Cerium atom is very sensitive to the 
crystalline and chemical environment and the differences 
in the interaction strengths of the localized 4f states with 
the conduction band can give rise to a wide variety of be- 
haviors. A lot of effort has been devoted to implement 
realistic calculations using different approaches P, [3, S| • 

In this work, we apply a mixed ab initio many-body 
technique to study the monopnictides CeX (X=P, As, 
Sb, Bi), the monochalcogenides CeY (Y=S, Se, Te) and 
CeXi_a;Ya; alloys which exhibit interesting and unusual 
physical properties. These compounds crystallize in the 
simple rock-salt structure making it possible to study in 
depth the influence of the electronic structure on their 
properties. In particular, the heavier compounds, namely 
CeSb, CeBi and CeTe 0] present a strong magnetic 
anisotropy, great sensitivity to the application of pres- 
sure and to the dilution of Cerium by non magnetic ions, 
as well as to the substitution of a pnictogen by another 
pnictogen or by a chalcogenjl|. These heavier compounds 
also show, in particular, anomalous small values for the 
4f crystal-fleld splittings (Acf) as compared to the other 
compounds within the series. We are interested here in 
doing a thorough analysis of the crystal field splittings 
as well as of the symmetry of the 4f states in these Ce 
compounds, for which is necessary to disentangle the in- 
teractions among the 4f states and the environment. We 
focus, thereafter, on the study of the crystal field splitting 
within the CeX and CeY series and also on its evolution 
along the CeXi-xYx alloys as a function of x. 

The anomalous value of the crystal field splitting of 



the heavier CeX and CeY compounds has been previ- 
ously treated and understood by Wills and Cooper 111, i^. 
These authors showed that the dominant contribution 
to the splittings in these monopnictides can be obtained 
from the point-charge (PC) model which is appropriate 
for insulators and ionic systems. But the depression of 
the splitting along the series can only be understood if hy- 
bridization effects are taken into account. The PC model 
accounts well for the splittings of rare earth monopnic- 
tides when the rare earth goes from Pr to Tb, but it fails 
to describe Ce systems because in these compounds the 
4f-band hybridization with the conduction states cannot 
be neglected. Wills and Cooper considered that the total 
splitting in these systems is the result of two indepen- 
dent contributions: the extrapolated value from the PC 
model considering non-hybridized 4f levels and the split- 
tings induced by hybridization. They calculated the 4f 
hybridization function out of the conduction band den- 
sity of states obtained using the linear mufiin-tin-orbital 
(LMTO) method within the atomic-sphere approxima- 
tion (ASA) for the self consistent potential. This last 
approximation does not consider the anisotropy of the 
crystalline environment. Within these calculations the 
authors treated the 4f state as belonging to the core . 
They obtained the hybridization contribution to Acf on 
the basis of the Anderson hamiltonian solved to second 
order in perturbation theory. 

In the present work, we obtain the Acf 's in a dif- 
ferent way, we perform ab initio calculations using the 
Full Potential-LAPW method within the local density 
approximation (LDA) [3| and treat the 4f states as part 
of the valence band. We then compute the 4f hybridiza- 
tion function following ref.0, Q- This function contains 
detailed information on the electronic structure of each 
system and is used as input for the Anderson impu- 
rity hamiltonian which is solved within a multiorbital 
non-crossing approximation (NCA). This LDA-NCA ap- 
proach has already been applied to cubic and tetrag- 
onal Ce based systems to calculate Ac^'s and to ob- 



2 



tain trends in the Kondo energy scales, yielding results 
which are in good agreement with the available experi- 
mental information [2, [§| . We calculate the crystal field 
splittings for CeX and CeY and compare them with the 
previous theoretical results and with experimental data. 
We also follow the evolution of A^f with concentration 
for the CeSbi_a;Te£!; alloys, which had not been previ- 
ously theoretically approached, and compare the results 
with experimental data. We also analyse the behavior 
of CeAsi-^^Se^;, which has not been previously tackled 
either experimentally nor theoretically. 

In section 2 we give a brief description of the LDA- 
NCA technique. Section 3 is divided into three parts: a) 
results for the Acf splittings of the CeX monopnictides, 
b) Comparison of the Acf splittings of the two series, 
CeX and CeY , and c) study of the evolution of Acf 
for the CeXi-xYx alloys. Finally, in section 4 we discuss 
and conclude. 



II. METHOD OF CALCULATION 



The i7 — > oo limit of the Anderson impurity model 
is solved by using the slave boson technique within the 
non crossing approximation (NCA). The NCA equa- 
tions consist of a couple of integral nonlinear equations 
for the pseudo-boson (/°) and pseudo-fermion (/^) self- 
energies, each of them containing the hybridization func- 
tion Tmm'i^)- This ?7 — > oo limit is reasonable to obtain 
crystal field splittings, because the splittings are stable 
within a wide temperature range above the Kondo tem- 
perature. For a detailed review of the NCA formalism 
see Ref 

With the LDA-NCA technique one can calculate the 
crystal-field splittings and also obtain the symmetry of 
the ground state and excited 4f levels. The crystal-field 
splittings are read from the separation of the peaks of the 
different spectral functions, Pmm's, which are shifted one 
with respect to the other due to the different degree of 
hybridization of each 4f level with the conduction band. 
We focus on the value of the splitting corresponding to 
the J = ^ multiplet, namely Acf = sjFy — Sfr^- 



The LDA-NCA approach is an ab initio many body 
technique that solves the Anderson impurity problem us- 
ing as input the hybridization function, r(e), of the con- 
duction band with the 4f-state of Ce. The hybridization 
function is calculated from first principles whithin the 
Density Functional Theory. In this work, the ab initio 
calculations are done using the full potential linearized 
augmented plane waves method (FP-LAPW), as imple- 
mented in the Wien2k code [lo| . 

For Ce systems the Anderson impurity hamiltonian has 
the form: 



H 



Lcka + ^£fmfLfm + U ^ Hfr. 



+ V {Vka,mfLcka + H.C.), 



(1) 



where Vka-,m is the hopping matrix element between the 
conduction electron states, (cka), and the 4f orbitals 
(/m), £fm is the corresponding 4f-energy level with re- 
spect to the Fermi energy and U the on-site Coulomb 
repulsion of the 4f electrons. r(£) is proportional to the 
product V^^ ^Vka,m' and, as suggested by Gunnarsson et 
al. 0, it can be estimated from the projected LDA 4f 
density matrix pf^^ at the Ce site in the following way. 
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In all cases the labels m and m' correspond to the 
different irreducible representations of the 4f states at 
the cubic Ce site. That is, for J = | the doublet Tr and 
the quartet Fg, while for J = | the doublets Fg and Fy 
and the quartet Fg. 



III. RESULTS 

The LDA calculations are performed at the experi- 
mental volumes for the CeX and CeY compounds. The 
muffin-tin radii, Rmt are taken equal to 2.4 a.u. in the 
case of Ce, while the corresponding radii for the anion 
ligand varies from 1.6 a.u. to 2.8 a.u. depending on the 
atomic radius. 102 k points in the irreducible Brillouin 
zone are considered to be enough for the quantities to be 
calculated. 

In the NCA-equations the 4f state has a bare energy 
value, which we take from photoemission experiments. It 
is namely, -3 eV for the series CeX , -2.6 eV for CeTe and 
-2.4 eV for, both, CeS and CeSe. All these energy levels 
are given with respect to the Fermi energy [l3| . 

Within NCA the self energies and the Green's func- 
tions are self-consistently obtained by considering the 
spectrum up to 6500 K above the Fermi level for the 
more delocalized systems (Y=S,Se, X=P,As) and up to 
5000 K for the more localized ones (Y=Te, X=Sb). 



A. Crystal-field splittings in CeX 

We discuss in this section the main features of the hy- 
bridization functions along the CeX series in order to un- 
derstand the evolution of the crystal-field splittings. As 
a representative example, we show in Figure [1] the calcu- 
lated hybridization functions for the Fy and Fg symme- 
tries corresponding to CeAs. In the inset, the detailed 
structure within an energy window around the Fermi 
level is given. 

It can be observed that the hybridization functions are 
very rich in structure. Far below and far above the Fermi 
level the Fg symmetry is the one with the largest hy- 
bridization contribution, while close to the Fermi level 
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the Ty one is stronger. This feature, with varying rel- 
ative weights, holds on for all the systems under study. 
The different relative weights are the finger prints which 
determine the evolution of the crystal-field splittings. As 
it will be shown in the next section, the 4/ states with 
r7 symmetry hybridize mainly with the Sd-Ce band while 
the ones with Fg symmetry do it mainly with the anion 
p states. 



reduction is more pronounced for the Fy than for the Fg 
symmetry. This decrease in the hybridization difference 
between both symmetries gives rise to a decrease in the 
value of the splitting. That is, the observed jump in the 
splittings along the series is correlated with a jump in 
volume when going from CeAs to CeSb. This analysis 
also applies to the CeY series (Ce monochalcogenides) . 




E(Ry) 

Figure 1: Hybridization function of the 4f states with the 
conduction band for CeAs : Fy (solid curve) and Fg (dashed 
curve) symmetries. In the inset, a zoom of an energy range 
around the Fermi level is given. 



B. Crystal field splittings in CeY vs CeX 



Due to the extra p-electron of the Y-anion, the Fg sym- 
metry is more strongly hybridized in the CeY series than 
along the CeX one. The Fy hybridization is of the same 
order of magnitude in both series, as it can be drawn 
from the comparison of Fig. [2] and Fig. [31 In these 
figures CeS and CeP are shown as examples. The same 
behavior applies for the other compounds. 

In the inset of figures [2] and [3] the partial p— and 5d— 
densities of states are plotted in detail. Comparing these 
densities of states with the corresponding hybridizations, 
it can be observed that below Ep, Fg hybridizes essen- 
tially with the p-states of the X and Y anions, while Tj 
does it with the 5d-states of the neighbouring Ce atoms. 



Due to the ionic character of the Ce compounds under 
study, the LDA unoccupied part of the spectrum is par- 
ticularly ill given. Taking this into account, in in order 
to reproduce the trends of the Ac_f's, we consider in this 
work the energy spectrum up to an energy of the order of 
6000 K above the Fermi level. Together with the correct 
trends in the splittings, we also obtain the experimental 
ground state symmetry, namely the Fy, for all these sys- 
tems. These results show that for these compounds the 
energy spectrum around the Fermi energy determines the 
symmetry of the ground state, even if the Fg hybridiza- 
tion function is the strongest one in the average. 

The calculated values for the splittings are: 165 K for 
CeP, 155 K for CeAs, 70 K for CeSb and 50 K for CeBi. 
The obtained values and trends are in agreement with 
experiments and with the previous theoretical work. The 
calculations reproduce properly the sharp decrease of the 
splittings when going from CeAs to CeSb. 

In order to understand the depression of Acf in CeSb 
and CeBi as compared to CeP and CeAsp^], we perform 
a calculation for CeP but at the experimental volume of 
CeBi, which is the largest one in the series. The effect 
of the negative pressure is, as expected, a reduction 
in the strength of the hybridization functions due to 
larger interatomic distances. Around the Fermi level this 
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Figure 2: Calculated hybridization functions for CeS: F? sym- 
metry (solid curve) and Fg symmetry (dotted curve). Inset: 
5d— DOS for Ce (solid curve) and p— DOS for the anion S 
(dotted curve). 
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Figure 3: Calculated hybridization functions for CeP: Ft sym- 
metry (solid curve) and Vg symmetry (dotted curve). Inset: 
5d— DOS for Ce (solid curve) and p— DOS for the anion P 
(dotted curve). 
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Table I: Experimental and LDA-NCA results for the crystal- 
field splittings, AcF=E(rs) - ^(Fy) , in the CeX and CeY 
series. The splittings are given in Kelvin. See Ref [H, [iBl for 
the experimental data. 

The values of the calculated Acf^s, shown in Table U 
are in very good agreement with the experimental results. 
The order of magnitude as well as the evolution of the 
Aci?'s are similar in both series. It should be noticed 
that the values corresponding to the CeY compounds are 
sHghtly smaller than the ones of the CeX series, and this 
is well reproduced in the calculations. 

In spite of the larger values of the Fg hybridization 
function for energies lying far below the Fermi level, the 
splittings are determined mainly by the hybridization 
character near the Fermi energy. The extra electron of 
the CeY compounds, as compared to the CeX ones, does 
not change the symmetry of the ground state nor that 
of the excited ones, preserving the order of magnitude of 
the Acf's when going from CeX and CeY. 

The feature which differentiates the two series under 
study, namely, the sHghtly smaller splitting observed in 
the CeY series as compared to the CeX one, can be corre- 
lated to the behavior of the partial 5d densities of states 
at the Fermi level. As it can be observed in the inset 
of Fig. [2] and [31 the Fermi level falls in a maximum of 
the 5d-D0S for CeP and in a valley of the 5d-T)0S in 
the case of CeS. The p-DOS, on the other hand, does not 



change considerably when going from CeP to CeS. The 
larger value of the 5d-D0S near Ep implies a larger Fy 
hybridization, which in turn gives rise to a larger split- 
ting for CeP than for CeS. This behavior of the partial 
5d-D0S at the Fermi level also holds for the other com- 
pounds. 



C. From CeX to CeY: Crystal-field splittings in 

the CeXi-iYx alloys 

So far, until now we have made a comparative analysis 
of the evolution of the Aci?'s along the two series CeX 
and CeY, the difference between them being the extra p- 
electron of the Y atoms as compared to the X ones. The 
Ac_f's of the solid solutions CeSbi_2;Te;r were measured 
by inelastic neutron scattering by Rossat [16,] for differ- 
ent Te-concentrations ( x= 0.05, 0.1, 0.5 and 0.7). The 
experimental data are shown in the upper plot of Figure 
m To the best of our knowledge, there is no previous 
theoretical work interpreting the shown results and com- 
plementing the suggestions done by Rossat in the sense 
that the observed trends in the alloys should be associ- 
ated to p — f mixing. 

In order to understand the non linear evolution of 
AcF as a function of Te concentration, we simulate the 
substitution of Sb by Te, and viceversa, by doing elec- 
tronic structure calculations within the virtual-crystal 
approximationfirl [Hi (VCA) in the two concentration 
limits. We have also analysed the splittings of the Hghter 
CeAsi_a:Sea: alloys within the same approximation. To 
this purpose, we add a small amount of valence electrons 
to the X atom of CeX and withdraw small amounts of 
valence electrons from the Y atom of CeY. Charge neu- 
trality is preserved in the calculations. Nearly one third 
of the added charge goes to the p levels of the anions, 
while the rest goes to the interstitial region. The occu- 
pation of the 5d and 4f levels of the Ce atoms, remains 
unchanged. Within this approach the original crystal 
symmetry is preserved. 

In the case of CeSbi-ajTe^;, we calculate the values 
of the splittings for four different VCA concentrations, 
namely by adding 0.1- and 0.2-electrons to Sb and by 
subtracting 0.1- and 0.2-electrons from Te in the par- 
ent compounds CeSb and CeTe, respectively. In or- 
der to compare with the behavior in the lighter alloys 
CeAsi_j;Sej;, we consider for these last systems two con- 
centrations, one adding 0.2e to As and the other sub- 
stracting this same amount of electrons from Se. In 
all cases the VCA calculations are done at the experi- 
mental volume of the corresponding parent compounds, 
as we are just studying here the electronic contribution 
to the evolution of the Acf values. In the first case, 
the obtained splittings are AcpiCeSb + O.le) = 90 K, 
AcriCeSb + 0.2e) = 120 K, AcriCeTe - O.le) = 55 K 
and AcpiCeTe — 0.2e) = 100 K. For the second simu- 
lated alloys the values are AcpiCeAs + 0.2e) = 200 K 
and AcF{CeSe — 0.2e) = 200 K. According to the exper- 
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Figure 4: Calculated Acf's of CeSbi-xTex systems. We have 
denoted CeSb+0.1 and CeSb+0.2 as a; = 0.1 and x = 0.2 and 
CeTe-0.2 and CeTe-0.1 as a; = 0.8 and x — 0.9, respectively. 
Experimental and calculated results are given in Kelvin. See 
Ref. [l^ for experimental data. 



imental trends for CeSbi-ajTea;, we obtain an increase 
in the value of the crystal-field splittings when a small 
amount of charge is added to CeSb and also when it is 
taken away from CeTe, as it is shown in Figure [H For 
the lighter alloys there is no available experimental re- 
sults for the splittings, but we obtain the same behavior 
as for the other ones. 

When CeX is doped with a small amount of charge the 
value of the Fg hybridization function increases slightly 
for energies below Ep. This is due to an increasing 4f- 
p mixing in that energy region. This extra charge also 
affects indirectly the 4f-5d hybridization through an in- 
creased p-5d mixing. This indirect mechanism induces a 
slight increment of the V^ hybridization function around 
Ep. There is no considerable contribution to the crystal- 
field splittings from energies lying more than 0.03 Ryd 
below Ep, because the hybridization functions of both 
symmetries change by a similar amount. This can be 
seen in Figure [H However, above Ep, the Fg hybridiza- 
tion goes down due to a decreasing number of unoccu- 
pied p states, while the Fy hybridization function goes 
up, these two effects give rise to an increase in the value 
of Aci?. These results agree with the experimental ones 
obtained by Rossat et al. and the output of our calcula- 
tions reinforce the interpretation done by them. 

In the other concentration limit, when a small amount 
of charge is removed from Te in CeTe, there is, in the av- 
erage, a slight decrease of the hybridization strengths be- 
low the Fermi level for both symmetries , namely Fy and 
Fg, but Fg is the one which shows the largest decrease. 
This last effect is to be attributed to the fact that Fg 
comes mainly from p-4f mixing and, that the withdrawn 
charge is essentially of p character. This growing dif- 
ference between both hybridization strengths is the rea- 
son for the experimentally observed larger values of the 
crystal field splittings, as compared to the splitting cor- 
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Figure 5: Frand Fg hybridization functions around the Fermi 
level for: CeSb {upper pannel) and CeSb+0.2 (lower pannel). 




Figure 6: Fyand Fg hybridization functions around the Fermi 
level for: CeTe {upper pannel) and CeTe-0.2 {lower pannel). 



responding to CeTe. In Figure [6] it can be seen that F7 
is larger than Fg below Ep for the system CeTe — 0.2e 
(lower pannel). In this plot, we indicate the symmetry 
separation (between F7 and Fg) with an arrow. 



IV. DISCUSSION AND CONCLUSIONS 

We apply in this work an ab initio-many body tech- 
nique (LDA-NCA) to calculate the crystal-field splittings 
of the monopnictides CeX and the monochalcogenides 
CeY series and their alloys. Within the abinitio frame, 
in which the 4f states are treated as part of the valence 
band, we calculate the hybridization function Fmm'(e), 
which contains detailed information on the electronic 
structure of each material. This Tmm'i^) function is used 
as input for the Anderson impurity model hamiltonian, 
which is solved within a multi orbital non-crossing ap- 
proximation in the U 00 limit. 

We show that the LDA-NCA technique gives the cor- 
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rect trends for the evolution of the crystal-field splittings 
of these systems, being the results in good agreement 
with the available experimental data and with previ- 
ous theoretical calculations We also obtain the cor- 
rect symmetry for the ground state of the 4f multiplet, 
that is, the Fy one for all the compounds and alloys. In 
particular, we correlate the sudden decrease of the ^cf 
when going from CeAs to CeSb (which are isoelectronic), 
to a big volume change. 

The most important contribution to the values of the 
crystal-field splittings comes from energies around the 
Fermi level. We obtain that the hybridization strength 
around Ep is due mainly to the 5d-4f mixing. In those 
systems where the 4f states are more locaHzed (CeBi, 
CeSb, CeTe) the intensity of the 5d-4f and of the p-4f 
hybridizations is smaller than in those where they are 
more delocalized. Due to this facts, magnitudes that 
depend on them, such as the crystal field splittings, are 
atenuated with respect to those compounds in which the 
4f states more delocalized. 

When comparing the results obtained for CeX with 
those for CeY, we obtain that the outcoming crystal-field 
splittings of the first series are slightly larger than those 
of the second one even if the CeY series has one more p 
electron. This is attributed to the fact that in the CeX 



series the Fermi level falls in a maximum of the 5d-D0S 
while in the CeY one it does it in a valley. The effect of 
the extra p electron in CeY is to produce an enlargement 
of the Fg hybridization far below the Fermi level, and this 
does not infiuence the value of the splitting considerably. 

Finally, we analyse the evolution of the splittings in 
the CeXi_a:Ya: alloys as a function of x. We consider 
one of the 'heavy' alloys, namely CeSbi-ajTe^;, and one 
of the 'light' ones, CeAsi_a;Se2;. Alloying gives rise to a 
non monotonic evolution of the Acf's which according 
to our results, obtained within the virtual-crystal ap- 
proximation, can be explained as an electronic effect by 
the only consideration of the added extra charge and to 
its effect on the hybridization functions. 
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